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ABSTRACT 





DNA barcoding is a powerful tool used for species identification by the expanding exploitation of mtDNA gene cytochrome c oxidase I (COI), as a genetic marker. In 
the present study, COI gene based DNA barcoding was employed to identify the epigeic species Eudrilus eugeniae. Molecular characterization of Eudrilus eugeniae 
was performed using Polymerase Chain Reaction (PCR) for COI gene. Phylogenetic tree was constructed using Molecular Evolutionary Genetics Analysis (MEGA7) 
software for the obtained sequences. The evolutionary divergences between its closely related species were also performed to disclose its amendments that occurred 
during evolution. 


INTRODUCTION: 

Research on taxonomy and biodiversity depends on efficient species diagnosis. 
DNA barcoding for animal species is based on exploitation of cytochrome c 
oxidase I (COI), which possesses high level of diversity. Mitochondrially based 
COI system of species identification provides much better taxonomic level of 
efficiency than morphological method of diagnosis [1]. DNA barcoding utilizes 
a short standard genomic portion referred to as 'DNA barcode'", which has suffi- 
cient sequence variation to discriminate among various species to achieve unam- 
biguous species level identification. The barcode may be constructed of one or 
several loci, together to be used as a complementary unit. One of the major 
advantages of DNA barcoding is that it provides 95% resolution of species-level 
recognition among mammals [2]. Also, it can be performed at any stage of life of 
an organism, whether live or dead [3]. 


Earthworms have the largest biomass (90%) among the soil invertebrates and 
constitutes a major component of the ecosystem. Disposal of sludge is an impor- 
tant yet tedious process in all industries. Among the cost-effective methods used 
for the treatment of sludge disposal, vermiculture is highly desirable. Eisenia 
foetida, Lampito mauritii and Eudrilus eugeniae together have been used in the 
biomanagement of papermill sludge [4]. The potentiality of these species was 
also assessed in terms of their efficiency and sustainability of vermicomposting 
water hyacinth, where Eudrilus eugeniae proved to be the most effective pro- 
ducer of vermicasts [5]. Identification of earthworm species is performed at an 
adult stage by dissection of the male genital organ [6], but that process is quite 
tedious, labor intensive and time consuming. Moreover, the other life stage 
worms remain unidentifiable with this method of identification. Hence, DNA 
barcoding is preferred as it applies to all life forms. 


MATERIALS AND METHODS: 

COLLECTION OF THE SAMPLE: 

Epigeic species Eudrilus eugeniae (Kinberg) was used in the present study for 
the degradation of the tannery sludge. Earthworm was collected from 
Ramancherry, Chennai Tamil nadu, India in plastic containers with soil and 
cowdung, transferred carefully to the laboratory with adequate moisture. They 
were handsorted and the species were identified on naked morphological obser- 
vations [7]. The biologically treated tannery sludge was collected from 
Pallavaram Tanners Association, Government Effluent treatment plant, Chennai 
Tamil Nadu, India. 


EXPERIMENTAL SETUP OF VERMI REACTORS: 

A Laboratory scale setup of vermi reactors was made up of plastic material in the 
form of open box provided with wire like holes at the bottom for drainage and ven- 
tilation. The reactors consisted of a basal layer of pebbles followed by coarse 
sand and a layer of garden soil and with a layer of dried cattle dung. These reac- 
tors were incubated with earthworms. The experimental units were maintained in 
triplicates and the control units were maintained five in numbers. In triplicates a 
total number of 36 experimental units were maintained to analyse the parameters 
at an interval of 12 days. The soil samples were collected from the triplicate 
experimental and control units of vermi composting reactors at an interval of 12 
days. The soil samples and earthworm samples were subjected to physio- 
chemical parameters DNA barcoding and sequencing. 


ISOLATION OF GENOMIC DNA: 

Isolation of genomic DNA was carried out by phenol chloroform extration 
method [8]. The pellet obtained in the last step was washed with 70% ethanol, air- 
dried completely, suspended in Tris-EDTA buffer and stored at -20°C until fur- 


ther use. The isolated DNA was assessed qualitatively and quantitatively by 
agarose gel electrophoresis and spectrophotometric method respectively. 


PCR AND DNA SEQUENCING: 

The polymerase chain reaction was carried out for amplification of the isolated 
DNA. The reaction volume was 20ul and the mixture consisted of forward and 
reverse primers, deoxyribonucleotide triphosphates (dATP, dCTP, dGTP and 
dTTP), template DNA sample and TaqDNA polymerase. Internal regions of COI 
gene were amplified using the following primers: COI-F(5'- 
GGTCAACAAATCATAAAGATATTGG-3') and COI-R(5'- 
TAAACTTCAGGGTGACCAAAAAATCA-3'). The reaction mixture was sub- 
jected to initial denaturation step at 94°C for 3minutes, followed by 35 cycles of 
denaturation at 94°C for 45 seconds, primer annealing at 47°C for Iminute, 
extension at 72°C for 1 minute and 20 seconds with the final extension at 72°C 
for 7 minutes. The PCR product underwent purification and DNA sequencing 
was performed by sanger sequencing (Applied biosytems 3500). 


SEQUENCE ANALYSIS: 

Nucleotide blast (BLASTn) using BLAST program and Genbank nucleotide 
database with default parameters was performed to determine the identity and 
the closest known relatives of the obtained sequences. Phylogenetic tree was con- 
structed using maximum likelihood method in MEGA version-5 (Molecular Evo- 
lutionary Genetics Analysis). The Distance Matrix Explorer, an action menu of 
MEGAS was used to compute the pair wise difference between the obtained tar- 
get sequence to its maximum aligned sequence. 


RESULTS: 

DNA was isolated from the earthworm sample and its purity was quantified by 
the A260/A280 ratio in a spectrophotometer. DNA concentration along with the 
purity details are given in Table 1. 


Table 1: Purity of the DNA observed in a spectrophotometer by calculat- 
ing A260/A280 ratio. 


Sample | Absorbance at) Absorbance at | Concentration 
260nm(A260) | 280nm(A280) (ng/ul) 


aiginy 
(A260/A280) 


Sample 2 0.128 0.070 6400 


Fig. 1: Genomic DNA isolation Eudrilus eugeniae. 
Lane 1: DNA from Eudrilus eugeniae, Lane 2: 1 Kb DNA Ladder 
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Fig. 2: COI gene amplification. 
Lane 1: COI from Eudrilus eugeniae, Lane 2: 1Kb DNA Ladder. 





Amplicons after PCR amplification were of 600bp for COI gene. Present study 
was conceded to determine the evolutionary relationship of Eudrilus eugeniae. 
The sequence obtained from the purified PCR product was subjected to compari- 


son with the nucleotide database BLASTn and was found to have maximum iden- 
tity to Vignysa teres, belonging to the family Hormogastridae. Construction of 
relationships between sequences is achieved using phylogenetic analysis that is 
aimed at identifying species and sequences from unknown samples [9]. Target, 
Eudrilus eugeniae branches to Vignysa teres, indicating the closest relationship, 
while it is distantly related to Megascolecidae species. 


KYBBE135 | Amynthas sp. JPA.2017 isolate APSAC E2 cytochrome cundase subunit | (COI) gene pattal cds mtachondnal 
HOMOSTTS. 1 Lumbnona sp. GL-2010 isolate T cytochrome ¢ cuodase subunt | (COM) gene pattal cds mtochondna 
——— KUBES IMS | Metaghwe sp EWEC.2 voucher E685 cytochrome ondase subunit 1 (00) gene pattal cds mitochondnal 
———— AB53231 1 Polygheretima tagrobanae mtochondnal CO! gene fr cytochrome onodase subunt | partal cds solate Ptap-2105 CO) 
GUD1208 | Acanthodrindae sp. DPEWST7ES voucher EWS} 
——— LOI 18734 | Pontodnivs itorals mtochondnal COI gene for cytochrome codase subunt | pattal cds aolate DW0G2) 
Tone 
J2US564 1 Vignysa teres voucher SG4 cytochrome cnadase sutunt | (C01) gene partial cds mtochondid 


47 cytochrome ondase subunt 1 (CON) gene partial cds mtochondaa 


KCOSBS 1 Amynthas moms: solate E27 cytochrome cundase sudunt | (COI) gene pattal cds mtochondna 
KP 214572 1 Megascolecidae sp. 23 MBL-2015 haplotype 23 cytochrome ondase subunit | gene partial cds mtochondnal 


Fig. 3: Molecular Evolutionary Genetics Analysis of target Eudrilus 
eugeniae with sequence producing significant alignment from database 
for COI region. 





Eudnius eugeniae voucher SGO U2 2015 cytochrome oxdase subunit ! (COl) gene. partial cds, mitochondnal 


Sequence (0: KAB32073.1 Length: 644 Number of Matches: 1 











Range 1: 75 to 636 Genbank Graphics 
Score Expect Identities Gaps Strand 
1120 bite(606) 0.0 610/612(90%) O/612(0%) Plus,/Plus 
Que i! i reaaa Ai TLiTaha " THMTITHITITT LETH ead eo 
Medtdtl any ‘an | Ree U ERAS aEaN PETEEELTTT TERE 
Se9ct 25 TATAAGGCTT CTTA itee A Tea eCAG CCOGGTECTTTTC TAGGAAGEGACCA Ba 
ery Te recA 
SURGGeen. i lil ORO GREGOR SOGRERE | | 
Sbict 85  ACTCTATAACACTATCGTTACAGCTCATGCTTTICTAATAATCITTTTCCTIGTAATGCE 144 
ery 2eP aL Ritalin Miiiiiiittl 
niitith TERELG Fd PERG Heed | | 
Seyct (145 sotitirsaT Too loce§riTTogsaas mC TCCCACT AS §TACTGGGASCOCCCGACAT 25. 
Query 1a] SASCATTCCCCCGSCTAASTATTITASGATTTTGATTATTAC TOCTTCACTSATTCTCTT 24e 
PEE EEE 
sejet 265 S26CeaTTClClCGel TAAST SATII AaGsA TIT TIGATIATTALCIC CT CatLtTaaertCicit | |6(26e 
eeeportiliiiinnitioniemitinntiiniTTTnine 
EET ELELTT i] | : PLT UTI Puuddail 
Sojct 265 AGTTTEGTCAGCTGCAGTTGAAAAGGGTGCAGGTACAGGA AACTOTTTACCOACCAET 324 
Query Jel TSG MTTITTITITTTTTTITTTT LITT ad ArT TTICTCTCTICATCT 360 
LURE SS a eT WORRORRCROORRRRe. 
Sbict 325 TGCAAGAARTCTTC CTCATECG COTCAGTAGACCTSGLCATTTTCTCTCTICATCT 364 
Query 361 MMII IIs 620 
i] Pe EUT Steel wan ryade 
Spjct 385 TSoCagococaTCaTcTatitiy lAgaecadtt: BACTTTATTA MACAGtAATTAATATAeS aas 
Query 421 ATGATCCGETCTTOGS aT Taenee niin neniian Ser TART ARTEMIS ae 
Ui | Peed eed bya 
SS3yct “45 TENN bit TICGATTAG4AaC katte ALTATTIGTATGSRGCAGT SOT SaTCaCTéT See 
query 461 SOTOCTRCTICTICTATCACTT CCARTCT TAG rf SGCAA hamn ACTICTCACAGA 54¢ 
—_ OEE SRGRaEE Hn 
Sect s85 ABISlraclTiTcrircrsrcaci ge CAGTCTTAGC LLabe ht ea Ape eCTCaCaGe 8 Sed 
Query 541 sp ainehthintiiall billinyainniih eli aanitittinien Bo 
BESERORE Pe 
S 5 "GAL TGEGGGTGGAGATCCAATICTITATICA 422 


CTTT 
eV He 
S7ct 565 FetAUt ead et tat rerteebere any 
T 


Fig. 4: 


Que Ol niin S 612 
1 | eae | 
Seict 625 ACATTTATICTS &36 
DISCUSSION: 


DNA barcodes have been proposed as a faster approach to provide novel species 
identification and discovery, unlike the taxonomists who would have spent 
decades to discover 10-15 million species using current description and publica- 
tions [10]. The highly conserved mitochondrial enzyme cytochrome oxidase c is 
coded by multiple genes containing regions that evolved at different rates [11]. 
These markers were exploited as DNA barcodes because of their potential to iden- 
tify putative regulatory elements, as they possess sufficient sequence diversity, 
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individually or in combination, to discriminate among species [12]. The usual 
and recommended method of DNA barcoding involves polymerase chain reac- 
tion (PCR) amplification of suitable regions of the genome, sequence analysis of 
the amplicons, and alignment of the analysed sequence with reference sequences 
[13]. 


Figure 3 exemplifies the Maximum likelihood tree constructed for the target 
Eudrilus eugeniae using MEGA7 software. The evolutionary history was 
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inferred by using the Maximum Likelihood method based on the Tamura-Nei 
model [14]. The tree with the highest log likelihood (-6600.24) is shown. The 
bootstrap consensus tree inferred from 500 replicates is taken to represent the evo- 
lutionary history of the taxa analyzed [15]. Branches corresponding to partitions 
reproduced in less than 50% bootstrap replicates are collapsed. The percentage 
of trees in which the associated taxa clustered together is shown next to the 
branches. Initial tree(s) for the heuristic search were obtained automatically by 
applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances 
estimated using the Maximum Composite Likelihood (MCL) approach, and then 
selecting the topology with superior log likelihood value. The tree is drawn to 
scale, with branch lengths measured in the number of substitutions per site. The 
analysis involved 10 nucleotide sequences. Codon positions included were 
1st+2nd+3rd+Noncoding. All positions containing gaps and missing data were 
eliminated. There were a total of 595 positions in the final dataset. Evolutionary 
analyses were conducted in MEGA7 [16]. The genomic analysis should that 
there were no major changes in nucleotide sequences after the earthworm was 
treated with tannery sludge for vermicompost. 


CONCLUSION: 

DNA barcoding has great potential for enhancing ecological and evolutionary 
investigations if the right genetic markers are exploited. In addition, phylogen- 
etic analysis is also important for identification of the potential changes in the 
nucleotides of a species and its variations. In the present study, molecular charac- 
terization of the earthworm Eudrilus eugeniae was carefully studied using COI 
gene and its evolutionary relationship was constructed using MEGA 7. 
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